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We study propagation of the Gaussian beam of spin waves and its reflection from the edge of thin 
yttrium-iron-garnet film with in-plane magnetization perpendicular to this edge. We have performed 
micromagnetic simulations supported by analytical calculations to investigate the influence of the 
surface magnetic anisotropy present at the film edge on the reflection, especially in the context of the 
Goos-Hanchen effect. We have shown the appearance of a negative lateral shift between reflected 
and incident spin wave beams’ spots. This shift is particularly sensitive to the surface magnetic 
anisotropy value and is a result of the Goos-Hanchen shift which is sensitive to the magnitude of the 
anisotropy and of the bending of the spin wave beam. We have demonstrated that the demagnetizing 
field provide graded increase of the refractive index for spin waves, which is responsible for the 
bending. 

PACS numbers: 75.30.Ds, 75.30.Gw, 75.70.Rf, 75.78.Cd 


In recent years magnetic nanostructures with con¬ 
trolled magnetization dynamics have been considered as 
candidates for design of new miniaturized devices with 
enhanced performance and functionality for various ap¬ 
plications, e.g. heat transport, energy conversion, mag¬ 
netic field sensing, information storage and processing.EH^l 
Spin waves (SWs), being propagating collective excita¬ 
tions of the magnetization are also regarded as informa¬ 
tion carriers, which can be exploited for information pro¬ 
cessing in devi ces p otentially competitive with standard 
CMOS systems.!^ Thus, understanding of SW properties 
in nanostructures is crucial in designing magnonic units 
and this is one of the main goal in the research field called 
magnonics.l^ It is expected that magnonic devices al¬ 
low energy-efficient processing of information which will 
combine the advantages of photonics (high frequency and 
wide band) and electronics (miniaturization) in a single 
unit.® One of the basic phenomena connected with wave 
propagation is the wave transmission and refiection.®tlil 
The reflection of SWs is determined by magnetic prop¬ 
erties of the film and boundary conditions at the border 
of the ferromagnetic material. The reflection of SWs has 
already b een in vestigated in theoretical and experimen¬ 
tal paper^l2*i2l where SWs were treated as plane waves. 
Use of wave beams, instead of the plane waves or spher¬ 
ical waves, in many cases, can be much more useful and 
opens new possibilities due to its coherence and low di¬ 
vergence. The known example of the wave beam is a 
light beam emitted by laser. Usually, its intensity pro¬ 
files can be described by Gaussian distribution (beams 
with such property are called Gaussian beams). How¬ 
ever, in magnonics the idea of SW beams is unexplored. 


with only a few theoretical and experimental studies con¬ 
sidering formation of SW beams at low frequencies due 
to the caustic or nonlinear effects.l^^^i^ 

An interesting phenomenon characteristic for the re¬ 
flection of beam is a possibility for occurrence of a lat¬ 
eral shift of the beam spot along the interface between 
the reflected and the incident beams - this phenomenon 
is called as the Goos-Haenchen (GH) effect. The GH 
effect was observed for electromagnetic waves acous¬ 
tic waves ,1^ electron^^ and neutron waves Also for 
SWs this topic was investigated theoretically for the re¬ 
flection of the exchange SWs (i.e., high frequency SWs 
with neglected dipole-dipole interactions) from the inter¬ 
face between two semi-infinite ferromagnetic films.^ It 
was shown that for the observation of the GH shift an 
interlayer exchange coupling between materials is cru¬ 
cial. Recently, we analyzed the GH shift at reflection of 
the SW’s beam from the edge of the magnetic metallic 
(Cobalt and Permalloy) and magnetic dielectric yttrium- 
iron-garnet (YIG, Y 3 Fe 50 i 2 ) films.l^We showed that the 
GH effect exists for dipole-exchange spin waves and can 
be observed experimentally. The magnetic properties at 
the film edge were shown to be crucial for a shift of the 
SW’s beam. 

In this paper we analyze the SW beam reflected from 
the edge of the thin ferromagnetic film. We focus our 
study on the magnetic properties of the film’s edge and its 
contribution to the shift of the SW beam. We show, that 
measurements of this shift can provide information about 
the local values of the surface magnetic anisotropy, and 
thus also about the local magnetic properties at the edges 
of the magnetic film. Our attention is concentrated on 
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detailed investigation of the SW reflection from the YIG 
film, a dielectric magnetic material highly suitable for 
magnonic applications due to its low SW damping, which 
is the smallest among all known magnetic materials!^ 
Recent experiments have shown possibility of fabrication 
of ve ry th in YIG films (with thicknesses down to tens 
of nrrPZEIl)^ which can be patterned on nanoscale and in 
which the SW dyna mics can be controlled with metallic 
capping layersPSEHl The magnetic properties of the film 
edge influence SW dynamics, their significance increases 
with decreasing size of device and will play important 
role in spintronic and magnonic nanoscale devicesP^^^^ 
However, edge properties at this scale are hardly accessi¬ 
ble to experimental techniques. In this paper, we propose 
a tool for the investigation of the magnetic properties at 
the edges of thin ferromagnetic film, which exploits shift 
of SW beams’ spot at the reflection. 

The micromagnetic simulations (MMS) and the ana¬ 
lytical model of the GH shift are described in section]^ 
Comparison of the results emerging from the analytical 
model with MMS, development of the model of SW bend¬ 
ing and the discussion of the results are presented in sec¬ 
tion]^ The paper is summarized in section III 


I. MODEL AND METHODS 
A. Model 

We consider a thin YIG film with the thickness, = 
5 nm, much smaller than lateral dimensions of the film 
{Lz <C Lx, Ly) as it is shown in figureThe film is mag¬ 
netically saturated by an in-plane static external mag¬ 
netic field H (we assume a value fXoH = 0.7 T) which is 
applied along the y-axis, perpendicular to the edge of the 
film. We study SWs which propagate in the film plane 
{x,y). The considered edge of the film is along x axis 
and located at y = 0. For description of the SW prop¬ 
agation, it is more convenient to define also the second 
coordinate system {x',y'). As it is shown in Fig. in 
this coordinate system the wave vector k, of the incident 
SW is parallel to y' axis and wave fronts are parallel to x' 
axis. Therefore, we can define the angle of incidence 9{ as 
the angle spanned between ki and normal to the edge (y 
axis). We limit angle of incidence to the value O-, = 60° in 
this study. We assume the SW frequency / = 35 GHz, at 
this frequency the propagation is almost isotropic in the 
film plane due to significant contribution of the exchange 
interactions (as confirmed latter in the paper with cal¬ 
culated isofrequency contours). In calculations we have 
used magnetic parameters for YIG at low temperatures: 
saturation magnetization Ms = 0.194 x 10® A/m and 
exchange constant A = 0.4 x 10“^^ J/m. An additional 
advantage of the YIG film is its relatively small static de¬ 
magnetizing field, which is proportional to Ms. All these 
properties of YIG simplify the analysis and helps us to 
focus mainly on the influence of the surface magnetic 
anisotropy on the reflection of SWs. The surface mag¬ 


netic anisotropy can have different origin, besides change 
of the crystallographic structure at the edge, the applying 
coating material and roughness can also influence surface 
anisotropy.^SlHlol Nonetheless, the microscopic mechanism 
of surface magnetic anisotropy is not the subject of this 
paper, our main concern is the influence of anisotropy on 
the reflection of a SW beam. 



Figure 1. Schematic plot of the thin YIG film geometry con¬ 
sidered in the paper. The film has thickness , which is much 
smaller than the film’s lateral sizes, Lx and Ly. The {x,y,z) 
coordinating system defines the structure with the film edge 
at y = 0 (hatched area). The coordinating system (x',y',z) 
defines the SW beam, with the wave vector parallel to y' and 
wave fronts parallel to x'. The area hatched by orange lines 
and located in the center of coordinate system {x',y',z) cor¬ 
responds to the excitation area, ki and kr are wavevectors of 
incident and reflected SW beams, respectively. AX is a total 
shift of the SW beam reflected at the edge. 


Magnetization dynamics is described by the Landau- 
Lifshitz (LL) equation of motion for the magnetization 
vector M: 


dM 

dt 


_ItL 

1 + 


iMxHoff — 


Mg (1 -|- a‘‘ 


.Mx(Mx H,ff), 
( 1 ) 


where: a - is the damping parameter, 7 - the gyromag- 
netic ratio, Hos - effective magnetic field. The first term 
in the LL equation describes precessional motion of the 
magnetization around the effective magnetic held and the 
second term enriches that precession by damping. The 
effective magnetic held in general can consist of many 
terms. In this paper we consider only the most impor¬ 
tant contributions: the external magnetic held H, the 
non-uniform exchange held Hex and the long-range dipo¬ 
lar held Hd: Heff = H -h Hex + Hd. 


B. Analytical model of the GH shift 

In our analytical study we consider SWs with large 
wavevectors, where the contribution from the dynamic 
dipole interactions is small and SW dynamics is mainly 
determined by the exchange interactions. In Eq. Q we 
will also make a linear approximation, which allows us 
to decompose the magnetization vector into a static part 
equal to the saturation magnetization Mg and a dynami- 
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cal part laying in the plane perpendicular to the direction 
of Ms: M = —Msy + ni{x,y,t). This approximation 
is valid when the dynamical part of the magnetization, 
|m| is much smaller than the saturation magnetization 
Ms- With this approximation we can assume harmonic 
time dependency for m oc where ui is the angular 
frequency of SW. SWs’ damping is neglected here. To 
study the incidence and reflection of SWs from the edge 
of the thin film, we start with dispersion relation for SWs 
in thin film assuming that the wave vector is in plane of 
the filmISl 

= [wH + + WM (1 - / (kL))] X 

[wh + ll^uJuk'^ + wm/ {kL) sin^ 6»k] , (2) 

where: y-o is permeability of vacuum, wh = 
is the internal field, for the in-plane magnetic field we 
assume Hi = H , wm = lyoMs, and the exchange length 
k is a wave number, 0k is the angle 
between the saturation magnetization and the in-plane 
wave vector k, and 


/ (a:) = 1 - 


X 


(3) 


We assume that the wave number is large, i.e., the con¬ 
dition khz 3> 1 is satisfied. Therefore, the SW disper¬ 
sion relation can be simplified (0k = in our geometry, 
see Fig. and reduced to the well-known Herring-Kittel 
equation: 


= (wH + ll^tOMk"^) (wH + ll^uJuk'^ + WM sin^ 0i) , (4) 


For further calculations we need formula for fc as a ex¬ 
plicit function of the frequency w. This dependence takes 
the following form: 


k^ 


AA 


-2H -Ms sin^ 0; -h 


I 4w2 


M| sin^ 


(5) 

The calculation of the reflection coefficient requires the 
introduction of the boundary conditions at the film edge 
(at y = 0) for m. We show below that the magnetization 
vector m components of the SW at the film edge (plane 
y = 0) satisfy the following boundary conditions of the 
Rado-Weertman typeP^ 


(y) 

dy 


dm (y) 


= 0 


y-O 


( 6 ) 


where d is an effective pinning parameter. This pinning 
parameter can take into account the dipole contribution 
in the finite width stripe {Ly is finite) and also contribu¬ 
tions from the exchange interaction and a surface mag¬ 
netic anisotropy at the edge of the film. The former con¬ 
tribution is anisotropic in restricted geometry and have 
to be calculated in our particular case. To do that we 
start from the boundary conditions written in the paper 
Ref. [23] (see Eq. (2) there): 


M X 


• dM 
dn 


-VM^-a + — 0 

do ) 


(7) 


is the energy density of the uniaxial surface anisotropy, n 
is the unit vector along the anisotropy axis, is a sur¬ 
face anisotropy constant, and Hj is the inhomogeneous 
magnetostatic field near the surface. The authors of the 
papers [431andl44] assumed that the static magnetization 
is parallel to the magnetic element surface and deduced 
the pinning of the dynamic magnetization component di¬ 
rected perpendicularly to the surface. Such perpendicu¬ 
lar magnetization component corresponds to increasing 
of the surface magnetostatic energy, which is an effec¬ 
tive easy plane surface anisotropy. An additional uni¬ 
axial surface anisotropy can be accounted leading to a 
re-normalization of the magnetostatic anisotropy. The 
system tries to reduce the surface magnetic charges de¬ 
veloping some inhomogeneous magnetization configura¬ 
tion near the surface. But this costs an additional vol¬ 
ume magnetostatic and exchange energy. The pinning 
was strong for thin (L^ = 2-20 nm) magnetic elements 
reflecting a balance between these energy contributions. 
The pinning calculated in Refs. |23| and |22] is a result of 
rapid change of the dynamic magnetostatic field near the 
surface (on distance about of L^). Whereas, in our case 
the static magnetization is perpendicular to the surface 
(film edge, the plane y = 0) and the dynamical mag¬ 
netization components are parallel to the surface plane 
xOz. Moreover, there is a SW wave vector component 
kx parallel to the surface. The pinning also has magne¬ 
tostatic contribution due to the strong dependence of the 
static dipolar field i4d (y) on the ^-coordinate near the 
surface plane y = 0. This dipolar pinning competes with 
the pinning induced by the surface anisotropy. There¬ 
fore, the contribution of even weak surface anisotropy is 
important in this case. We can re-write the boundary 
condition equations in the symmetric explicit form: 


iLdm 

Lz dn 


-f h-h 


1 

2 


Ks 


m = 0, 


( 8 ) 


where Kg = 2Ks/y,oMg, Hd = HJ) -1- h and the contribu¬ 
tions from exchange, dynamical dipolar, static dipolar, 
and surface anisotropy fields, respectively, are included. 
The magnetization precession in such non-ellipsoidal el¬ 
ement under consideration (see Fig. is elliptical due 
to non-equivalent dynamical dipolar field components hx 
and hz along the film normal Oz and in-plane Ox direc¬ 
tions. 

The dipolar field components can be expressed via the 
dynamical magnetization components using the method 
of the tensorial magnetostatic Green functions, see Ap¬ 
pendix 1^ and Ref. |2S| . Even assuming that the dy¬ 
namical magnetization does not depend on the thick¬ 
ness coordinate z and averaging over z we still have two- 
dimensional problem because the magnetization depends 
on the in-plane coordinates x^y. We assume that the 
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dynamical magnetization can be represented in the form 
ma{x,y) = exp {ikxx) ma (y), where kx is almost con¬ 
tinuous variable due to large element size along Ox di¬ 
rection. The magnetization profile TOq, (y) can deviate 
from the plane wave near the element edge, where dipo¬ 
lar helds are strongly inhomogeneous. Using the Green 
functions formalism we simplify the boundary conditions 
given by Eq. (|^ and write them in the generalized Rado- 
Weertman form: 


dnia 

dy 


+ daTTla 


= 0, where a = x,y 


y=o 


(9) 


but the pinning parameters dx and dy are different and 
depend on the wave vector component kx- We get 


dx 

dz 


Lz 1 / {kxLz') 

W. [ 2 Tz\ ’ 

Lz f f{kxLz) Ks \ 

IL I 2 LzJ- 


( 10 ) 


The typical value of kxL^ is order of 1 that corresponds 
to wave number kx of about 0.1 nm“^, i.e. to the dipolar- 
exchange SW regime. To simplify the further analytical 
consideration we assume that the exchange energy dom¬ 
inates and neglect the dynamical dipolar fields in the 
boundary conditions. Therefore the symmetric pinning 
parameters d = dx = dz are: 


d = 


P 

'■ex 



( 11 ) 


Due to translational symmetry along the interface {Ox 
direction), the wave vector component along the interface 
{kx) should be conserved. As a consequence, the angle 
of incidence is equal to the angle of reflection, 0, = 0,.. 
Therefore, based on Eqs. and (|^ and assuming ma {y) 
dependence in the plane wave form we can derive the 
equation for the reflection coefficient R (see Appendix[B|): 


i^/k"^ — kx + d 
\sjkp — kx — d 


( 12 ) 


When the incident SW is represented as a wave packet 
of a Gaussian shape, with a characteristic length in mo¬ 
mentum space Afc(. k', then according with the sta¬ 
tionary phase methocpiliSI the reflected beam reveals a 
space shift relatively to the incident wave packet of the 
length (the GH shift): 


AAqh = — 


dij) 
dk' ’ 


(13) 


where = arctan (3(i?)/3fi(ii)) is the phase difference 
between the reflected and incident waves, 3(i?) and 3?(i?) 
are the imaginary and real parts of the reflection coeffi¬ 
cient calculated from Eq. (121. Thus, the GH shift AATgh 
dependence on the magnetization pinning coefficient can 
be expressed by the following equation: 


AXqh 


2d tan 6i 
d? + {kcosOif 


(14) 


C. Micromagnetic simulations 


Micromagnetic simulations (MMS) have been proved 
to be an efficient tool for the calculation of SW dy¬ 
namics in various geometriesWe have exploited 
an interface with the GPU-accelerated MMS program 
MuMaxli^ which uses finite difference method to solve 
time-dependent LL Eq. ([T). 

In our MMS we consider SWs propagation in thin- 
fllms and reflection from the him edge. Simulations 
were performed for the system shown in Fig. [l] of size 
4000 X 12000 X 5 nm {Lx x LyX L^), which was discretized 
with cuboid elements of dimensions 2.5 x 2.5 x 5 nm 
{lx X ly X lz are much less than 13 nm, i.e., the exchange 
length of YIG). The surface magnetic anisotropy was in¬ 
troduced in MMS by uniaxial magnetic anisotropy value 
ATu in the single row of discretized cuboids at the him 
edge according with the dehnition = ATg/Zy,!^ where 
ly = 2.5 nm is a size of cuboid along normal to the edge. 

Simulations consist of two parts according to the al¬ 
gorithm presented in Fig. First, we obtain the equi¬ 
librium static magnetic conhguration of simulated sys¬ 
tem. In this part of simulations we start from ran¬ 
dom magnetic conhguration in presence of high damping 
{a = 0.5). Then, the results of the hrst stage are used 
in the dynamic part of simulations during which a SW 
beam is continuously generated and propagates through 
the him with reduced, hnite value of damping parame¬ 
ter, a = 0.0005, comparable with the values present in 
high quality YIG hlms.^^ The SWs are excited in the 
form of a Gaussian beam. After sufficiently long time, 
when incident and rehected beam are clearly visible and 
not changing qualitatively in time the data necessary for 
further analysis (‘POSTPROCESSING’) are stored. 



Figure 2. Algorithm of the SW dynamics calculations by 
MMS. MMS consist of two steps, in the hrst step the sys¬ 
tem is stabilizing-the equilibrium magnetic conhguration is 
obtained. In the second step using the stabilized magnetic 
conhguration SWs are generated by applying small rf mag¬ 
netic held. The data stored during MMS are processed during 
the stage called ‘POSTPROCESSING’-the hnal results of the 
SW dynamics are extracted. 

To generate Gaussian beam of SWs we introduce a nar¬ 
row rectangular area (excitation area, marked in Fig. 
with orange dashed lines) with a long side parallel to the 
expected wave fronts (along x' axis). Within the exci¬ 
tation area we introduce a radio-frequency (rf) magnetic 
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field oscillating at frequency of 35 GHz, h^yn {x',y',t) = 
/io(x', y') exp (iwt)). The field hdyn is perpendicular 
to the static magnetic field and its amplitude changes 
along the x'-axis according with the Gauss distribution 


G{x) = exp 



Xo'f 


and along y'-axis takes 


uniform, non-zero values only for small window of width 
w — 5 nm centered around yo' . I is the length of the 
excitation area (in our simulations I = 1.5 pm) and 
(a = 0.2) is a parameter which can be treated as variance 
of the Gauss distribution centered around xq'- Hence, 
ho {x', y’) = /i0H(y'-yo'+w^/2)0H(-y'+yo+w/2)G(a;'), 
where 0 h is the Heaviside step function. We assume that 
h = Q.Q2H being maximum amplitude of the rf magnetic 
field which needs to be small to stay in linear regime. 
Example result of MMS is shown in Fig. 


Extracting the precise value of the shift of SW beam 
AX from MMS results requires the following three-stage 
procedure. In the first stage we extract a series of SW 
intensity profiles using an array of screen detectors paral¬ 
lel to the y-axis for different locations (xj) along the film 
edge but far enough from the reflection point, i.e., out of 
the interference pattern of SW near the reflection point. 
At every point Xj the intensity was calculated using equa¬ 
tion: (y) = \m^ {xj,y,t)\dt, where T = 1/f and 

rUz is the component of the magnetization vector per¬ 
pendicular to the film plane. In the next stage, using 
Gaussian fitting, we have extracted positions of centers 
of the intensity profiles for every yj{xj). Having series of 
the peak positions and its locations along y (red and blue 
full dots in Fig. |^for the incident and reflected beams, 
respectively) we can extract rays of the incident and re¬ 
flected SW beams (red and blue solid line in Fig. re¬ 
spectively). Finally, the value of the shift AX can be 
easily calculated with small errors up to several nanome¬ 
ters. 


II. RESULTS AND DISCUSSION 


Result of the pure GH shift in dependence on the 
pinning parameter obtained from the analytical model 
[Eq. (14)] is presented in Fig. |^a). This dependence 
AXGH(d) is an antisymmetrical function with respect to 
d = 0 and has maximum and minimum value for d < 0 
and d > 0, respectively. The effective pinning param¬ 
eter d has contributions from the exchange interaction, 
dipole interaction and magnetic surface anisotropy. To 
study influence of the magnetic anisotropy we show also 
AX in dependence on Kg in Fig. ib) with solid line. 
GH shift exists for Kg = 0 (marked with small square in 
Fig.|^b)) due to effective pinning coming from the dipole 
interaction^^ and is AAgh = —17.0 nm. AAgh takes 
maximum absolute values for Kg = 0.271 mJ/m^ and 
Kg = —0.153 mJ/m^. AAgh = 0 for Kg = 0.59 mJ/m^, 
this is when magnetic surface anisotropy compensate the 
effect of the dipole interactions at the film edge. For large 
negative and positive values of Kg the GH shift tends 
monotonously to zero. This shows that the measure of 
the GH shift can be used to indicate the surface mag¬ 
netic anisotropy at the thin film edge locally, especially 
in the range of its sudden change, i.e., between —0.153 
and 0.271 mJ/m^. To test this possibility we perform 
MMS according with the procedure described in section 

ESI 



d [1/nm] 




Figure 3. Exemplary result of MMS showing the color map 
of the dynamic magnetization amplitude of the SW beam at 
reflection from the edge of the YIG film. The presented result 
was achieved for Kq = —0.1 mj/m^. At this value of the 
surface anisotropy the negative value of SW beam shift is 
observed: AX = —45.8 nm. The blue line corresponds to the 
rays of the incident beam and red line to the reflected beam. 
The dashed white line close to the ray of the incident beam 
points at the direction of the SW wavevector. 


Figure 4. Analytical results of the GH shift in the reflection of 
the SW beam from edge of the thin YIG film calculated using 
Eq. (141 in dependence on (a) pinning parameter d and (b) 
magnetic surface anisotropy constant Kg. The blue squares 
corresponds to values of GH shift for d = 0 in (a) and Kg = 0 
in (b). 


Dependence of the SW beam shift on the surface 
anisotropy constant obtained from MMS for yo77 = 0.7 T 
in thin YIG film is presented in Fig. with green solid 
dots. The value of the shift for Kg = 0 obtained 
from MMS is AA = —32.4 nm and this is signifi¬ 
cantly larger than the GH shift obtained from analyt¬ 
ical solutions (AAgh = —9 nm). The maximal value 
of AA = 13.83 nm is found for Kg = 0.25 mJ/m^ and 
minimal is AA = —49.6 nm for Kg = —0.2 mJ/m^, i.e., 
out of the scale presented in Fig. This dependence is 
quantitatively similar to the function obtained in the an¬ 
alytical model (Fig. 1^. Nevertheless, there are distinct 
differences between both results. 

The analytical model is based on number of assump- 
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Figure 5. Results of the MMS (green solid points) and 
obtained from the analytical models presenting dependence 
of the SW beam shift on the surface magnetic anisotropy 
constant in the YIG film in the external field 0.7 T. Only 
the results for fully saturated sample Ka > Kq {Ko = 
—0.22 mJ/m^) are shown. Dashed blue line corresponds to 
basic analytical model for GH shift [Eq. ( |14| )] without in¬ 
cluded SWs bending corrections (c = 1). Dash-dotted black 
line presents results for the analytical model [Eq. (201] with 
c = 0. Solid orange line correspond to results of the analytical 
model with fitted c = 0.82. 


tions which are absent in MMS, thus there is a dis¬ 
agreement between both results. The full saturation of 
the magnetization is important assumption made in the 
analytical model. Surface magnetic anisotropy at the 
edge of the film increases (or decreases) a value of a 
total internal magnetic field at this edge. This change 
of the internal magnetic field in the single computa¬ 
tional cell at the edge is equal to the surface anisotropy 
field ^qHs = 2Ks/ {MJy). Therefore, exactly at the 
film edge the demagnetizing field can be compensated 
or enhanced by surface anisotropy field, in dependence 
on the sign of Kg. When surface anisotropy constant 
Kg = Ko = -\y.oL^Mg{H - Mg) w -0.22 mJ/m^ the 
total internal field at the surface is equal 0 (the demagne¬ 
tizing field at the edge is compensated by the anisotropy 
field). Hence, for Kg < Kq the equilibrium orientation 
of the magnetization at the edge of the film rotates from 
the saturation direction (this is transformation from the 
easy axis to the easy plane configuration). This is demon¬ 
strated in Fig. §(b), where magnetization components 
along X, y and z axis are presented as a function of the 
distance from the film edge for Kg = —2.0 mJ/m^. We 
can see the rotation of the magnetization from the y di¬ 
rection (saturation direction) towards the x axis. The 
change of the magnetization configuration is not taken 
into account in the model developed in Sec. m More¬ 
over, in real sample there is a possibility for appearing 
domain walls along the film edge, which can create ad¬ 
ditional factor for complexity of the problem. Therefore, 
in this paper we limit the analysis to the fully saturated 


sample, i.e., when Kg > Kq (we note, that the exact 
value of Kq depends on the magnitude of the external 
magnetic field). 



Figure 6. (a) Effective magnetic field along the y-axis inside 
YIG film for Ka = 0 is shown with blue solid line, orange 
dashed line corresponds to the external magnetic field (0.7 
T). (b) Static magnetic configuration in the vicinity of thin 
YIG film edge for Ka = —2 mj/m^ (i.e., Ka < Kq); blue solid, 
orange dashed and green dash-dotted line marks x, y and z 
component of the magnetization vector normalized to unity, 
respectively. 


The model shown in Sec. HD doesn’t take into ac¬ 
count influence of inhomogeneity of the internal mag¬ 
netic field on the behavior of propagating beam espe¬ 
cially in the vicinity of the film edge. It is probably the 
next most important factor, which makes a comparison 
of the analytical and MMS results difficult for the sat¬ 
urate state. In MMS, and in a real sample, the magne¬ 
tization which is perpendicular to the film edge creates 
an inhomogeneous static demagnetizing field, which is 
directed opposite to the magnetization saturation. In 
result, the internal magnetic field decreases monotoni- 
cally when moving from the film center towards its edge 
[solid-blue line in Fig. [^a)]. In our interpretation this 
inhomogeneity is responsible for a shift of AX in MMS 
towards negative values as compared to the results of 
Eq. ( [T4| ) [see. Fig. [^: for very high positive Kg the 
AXqh monotonously tends to 0 in the analytical model 
[Fig.|^b)], while results of the MMS show that the value 
of AX reaches non-zero negative value even for very high 
value of Kg (AX = —15.2 nm for Kg = 10 mJ/m^). This 
means, that the inhomogeneity of the internal magnetic 
field in the close vicinity of the thin fil m edg e causes an 
increase of the refractive index for SWa^^EI] and conse¬ 
quently results in bending of the SW beam and changes 
the shift measured in far field. 

Here, we propose the simple analytical model of the 
wave bending, which allows to estimate factor which cor¬ 
rect the value of the GH SW beam’s shift obtained from 
Eq. (14). We will consider gradual change of the refrac¬ 
tive index for SWs in the vicinity of the film edge. 

Similarly as in optics, the refraction law for SWs 
(i.e., Snell law) can be concluded from analysis of the 
isofrequency contours and momentum conservation of the 
wavevector component parallel to the edge {k^ = const.). 
Let us assume, that full distance of the gradual change 
of the refractive index in thin film can be divided into 
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N thin slices, numbered with integer n. Therefore, 
we can assume multiple refractions on (N — 1) paral¬ 
lel planes, separating neighbor slices [Fig. [^b)]. At the 
plane between two arbitrary slices n and n -|- 1 accord¬ 
ing to Snell low [kn sin = fc„+isin6>i,„+i) [Fig. [3 a)]. 
Thus, we can calculate final angle of incidence after 
passing N slices: sindj^jv = ^sin0i_o, where and 
are wavenumber and incident angle in the last slice, 
/co = fc(yo) and 01 ^ are wavenumber and incident angle in 
the initial media, i.e., in the interior far from the edge of 
the film. Rewriting the expression for space coordinates 
(substitute k{y) = k^ and 9i{y) = 0i,N) we obtain: 


0i{y) = arcsin 



(15) 


From here, the final value of 9i{y) at the y distance from 
the him edge can be calculated, if the value of k{y), initial 
values of the wavevector fcg and the initial angle of the 
incidence are known. 


Assuming inhnitesimally small distance between inter¬ 
faces Ay —>■ dy we can transform summation into inte¬ 
gration along the y axis between initial point of the SW 
generation yo and the point of the rehection yjy: 


x{yo,yN) 



dy tan arcsin 




(17) 


Therefore, if the edge where reflection takes place is lo¬ 
cated at yAT = 0 (like in Fig. [^, the SW beam shift 
introduced by bending and measured in far held (obser¬ 
vation point is assumed at the same distance from the 
edge as the source of the SW) is described by following 
formula: 


— 2 [yo tanx(yo,0)], (1^) 




where factor 2 is included to take into account the beam 
way from source to the edge and after rehection to the 
observation point. 

The total shift of the SW beam observed in far held 
is a sum of the GH shift and the shift resulting from the 
bending: 

AA = AAbending + AAgH- (19) 


Figure 7. (a) Explanation of the refraction of the wave on the 
interface between media with low refractive index (n-th slice, 
orange line) and high refractive index (n -f 1 slice, blue line) 
based on isofrequency contours analysis. The conservation 
of the X components of the wavevectors in the refraction is 
required by the translational symmetry along the x. (b) Ex¬ 
ample of two refractions on the interfaces between n-th and 
n -I- 1, and between n -f 1 and n -|- 2 slice. The refractive index 
increases with increasing index n. The beam shift AAbending 
resulting from the bending is shown schematically. 


The knowledge of the incident angle at the vicinity 
of the film edge 9i{y) allows us to derive formula which 
describes propagation of the beam through the area of 
gradually changed refractive index. Similarly to previous 
approach, let us consider single refraction on the interface 
between n and n -I- 1 slice, as it is shown in Fig. Bb). 
Projection of the beam path onto the direction tangential 
to the interface is Axn+i = Ay tan0i_„_|_i where Ay is 
the thickness of the slice (we assume here for simplicity 
slices of the same thickness). Generalizing this formula 
for cascading refraction on N slices we obtain: 


XN 


N N 

Axn = Ay tan 

n—1 n—1 


N 

AyJ] tan arcsin 

n—1 



(16) 


Eq. (19) is a general formula describing a total shift of the 
wave beam propagating in media with gradual change of 
the refractive index. Way of the beam ray depends on 
relation describing k(y). And this is main bottleneck in 
this approach: unknown formula for wavevector of the 
SW in dependence on y in an area of the inhomogeneous 
effective magnetic field. However, knowledge of the dis¬ 
persion relation in homogeneous film, Eq. (^, can help 
in qualitative modeling of the total shift of the SW beam 
also in a part of the film with gradual change of the re¬ 
fractive index. 


Demagnetization field decreases a value of the in¬ 
ternal magnetic field, its dependence on distance from 
the edge can be expressed by the relation: H^{y) = 
4Ms arctan [L^/(2y)].l^ This field far from the edge of 
the film tends to zero [Fig. |^(a)]. We can substitute 
bias magnetic field in Eq. with the internal field 
H —)• H{y) = H — H^{y). This approach should be 
valid for slow change of i7d(y) with distance. However, 
in the vicinity of the thin film edge we observe a rapid 
change of the demagnetizing field [Fig. 6(a)] and the 
plane-wave approximation fails. Therefore, we propose 
to introduce homotopic transformation of the demagne¬ 
tizing field: H{y) ^ H - [cHd{yo) -b (1 - c) Hd{y)] with 
parameter c S [0,1], which will reduce influence of the 
rapid changes of the demagnetizing field on wave vector. 
Then, the wave vector magnitude can be described by 




















the following equation: 


kHy,ei) = 


ypMs ^ 


AA 


-2{H-[cHdiyo) + il-c)Hdiy)]) 


—Mg sin^ Oi 


+ M| sin^ 9i ), (20) 


where yo can be interpreted as position of the source of 
the SW beam, i.e., the position where the demagnetizing 
field magnitude is close to zero. 

Now, the only issue is a correct choice of the parameter 
c in Eq. (20) to find k{y), and to fit AX{Ks) dependence 


to the curve obtained from MMS. Taking c = 1 we as¬ 
sume, that isofrequency contours are the same in whole 
sample (fc does not depend on y), the wave propagate 
in homogeneous internal magnetic field equal H [because 
Hd{yo) — 0] and Eq. (20) reduces to Eq. ([^. In this 
case AJfbending = 0 and AX = AAgh, i-e, the total 
shift is equal to Eq. (14). The obtained curve is plot¬ 
ted in Eig. with blue dashed line. In opposite limit, 
c = 0 the dependence k{y) follows exactly the change 
of Hd{y). In this case the AAtotai(.?i^s) calculated from 
Eq. (19) is shown in Fig. with black dotted line. It 
takes value AX = A Abending = — 79 nm for large Kg, 
much below the value obtained from MMS. This discrep¬ 
ancy exists, because for each y we took in calculation of 
the AAbending [ln the integral Eq. 0] the dispersion 


relation Eq. (20) which is for the film with homogenous 


magnetic field. In real situation, the dispersion relation 
in the area of inhomogeneous refractive index will be dif¬ 
ferent from the local value. Thus averaging of the de¬ 
magnetizing field across some distance shall improve the 
estimation. Moreover, the c shall depend on the rela¬ 
tive value of the wavelength to the special changes of 
the refractive index. Thus, the value of c from Eq. (20) 


needs to be treated as a parameter, which includes an 
effective influence of the inhomogeneity of the refractive 
index on the dispersion relation of SWs. For c = 0.82 
we have obtained very good agreement between results 
of MMS and analytical model AX (Kg) « AX {Kg) for 
Kg > —0.22 mJ/m^, as it is shown in Fig. with orange 
solid line. 

The value of AX is very sensitive for small change of Kg 
between the extremes. The magnetic surface anisotropy 
in ferromagnetic films can take different values, how¬ 
ever the most interesting is the range around 0, where 
the transition from easy axis into easy-plane anisotropy 
takes place. Thus, the measure of the SW beam’s shift 
can indicate the local surface magnetic anisotropy at 
the film edge with spatial resolution limited by the size 
of the width of the SW beam. This information shall 
be important also for understanding and exploiting SW 
excitations and actuation at the surface of YIG film 
being in contact with Pt, where the magnetic surface 
anisotropy w as shown to play a significant role in the 
spin pumping .1^^^211111111 Further investigation is required 
to test an influence of the second ferromagnetic material 
attached to the YIG film edge on the SW reflection. Here, 


extension of the analytical model wi th pr operly defined 
boundary conditions will be required.^^IEIl 

The micromagnetic simulations were conducted for the 
value of the damping parameter a = 0.0005, which is 
close to the value of thiclP^ and also very thin films 
(tens of nm thick) of YIG as demonstrated experimen¬ 
tally recently.^ We have performed additional simula¬ 
tions for slightly smaller and larger damping to check the 
influence of damping on GH shift. Apart from increase 
of the amplitude of the reflected SW beam with decreas¬ 
ing a the results of the GH shift are very close to those 
presented in the manuscript. This indicates that small 
changes of the homogeneous damping does not affect GH 
shift. However, for the electromagnetic waves reflected 
from the interface GH shift is strongly affected if t he re- 
flected material is a medium with strong absorption j^^ l ^^ l 
where the absorption can even change the sign of the GH 
shift. Thus, we can suppose, that also in magnonics the 
inhomogeneous damping, especially with its high value at 
the interface or in the media behind the reflection edge, 
will influence GH shift. Further investigations are neces¬ 
sary to elucidate the role of inhomogeneous damping on 
the GH shift in the reflection of SWs from the edge or 
interface ferromagnetic films. 


III. CONCLUSIONS 


We have performed analytical and numerical study of 
the SWs beams shift at the reflection from the edge of 
the YIG thin film in dependence on the surface magnetic 
anisotropy present at the film edge. The GH effect and 
SW bending are shown to contribute to the SW beam’s 
shift measured in far field. We have shown that the GH 
shift is modulated in a broad range by changes of the sur¬ 
face magnetic anisotropy constant between two extremes: 
ATs^min = —0.155 mJ/m^ and ATg^max = 0.275 mJ/m^. It 
means that even small change of the surface magnetic 
anisotropy, some imperfections or changes in a surround¬ 
ing of the film edge can result in significant change in 
the reflected SW beams position in far field. The demag¬ 
netizing field gradually changes the refractive index of 
SWs of the film with its sudden increase near the film’s 
edge. This variation of the refractive index bends the SW 
propagating towards (and outwards) of the edge and in¬ 
troduce additional shift of the SW beam. However, this 
shift is independent on the surface magnetic anisotropy. 
For large positive values of Kg the value of AX is almost 
insensitive to changes of the magnitude of the magnetic 
surface anisotropy and is mainly a result of the SW bend¬ 
ing due to gradually increased refractive index at the film 
edge. These results shall be of importance for magnonics, 
its applications for sensing and also for developing a new 
direction of research devoted to metamaterial properties 
for SWs, especially the graded index magnonics. 
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Appendix A: Green functions 


The X and z components of the magnetostatic fields 
can be expressed as: 


r] = 0 because of the condition ^ 1- This justify us¬ 
ing the Taylor series decomposition near the point t] = 0 
calculating the dynamical dipolar fields and derivation of 
the boundary conditions given by Eqs. (101 and (11). 


hx{x,y) 

hz{x,y) 


dx'dy'GXX(x, y,x',y')mx{x’,y’), 
dx'dy'G^zix, y; x', y')mz{x', y'), (Al) 


Appendix B: Reflection coefficient 

The plane wave solutions are assumed for the dynamic 
components of the magnetization vector mQ.(r,t): 


where 


Gc.p{p-,p') = J^j d^\^Go^p{\^y^ G-P')^ (A2) 

G..(k) = -f{kL,)^, Gzz(k) = - [1 - f{kL,)], 


and function / is defined in Eq. (§. 

Substituting ma{x,y) = (ikxx)ma(]j) to Eq. (Al) 
and performing integration over x' and fc^-we can get the 
dipolar field in the form 


ha{x, y) = J dy'gaaiy, y')'ma{y'), (A3) 


where the magnetostatic kernels are calculated using 
Ref. [12]: 




by substituting this solution to the boundary condition 
([^ we obtain: 


+ d -k = 0. (B2) 

V /J v=o 


At the reflection point k^ ^ = k^,x = kx and ki y = 
—kr^y = ky we get: 


ikyAa — ikyBg + d{Aa + B^) = 0. (B3) 

Finally, the formula (12) for refractive index can be easily 
derived: 


Ba iky + d i^k'^ — k^ + d 

Aa iky — d i^kP' — k'^ — d 


(B4) 
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and 


-Xp) =-\kxQ at » 1, 


(A5) 


above, Kq is the modified Bessel’s function of the zero 
order, and y = y — y'. The zz-component of the mag¬ 
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